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Abstract 

o 

I The dromion of the Davey-Stewartson-1 equation is studied under 

' perturbation on the large time. 

a^ : 

^ I In this work we construct the asymptotic solution of the Davey-Stewart- 

O-f son-1 equations (DS-1): 

^' 1 

^ ■ tdtQ + -(<9| + + (Gi + G2)Q = etF, 

diG, = -^d,\Q\\ dr,G2 = ~d^\Q\'. 

■ Here e is small positive parameter, F is operator of the perturbation, the 

values of the parameter a = ±1 correspond so called focusing or defocusing 
DS-1 equations. 

The equations (|l|) at e = describe the interaction of long and short 
waves on the liquid surface if the capillary effects and potential flow are 
taken into account |1|, |^ . The theorems about the existence of the solutions 
for this equations in the different functional classes are known [Q, ^. The 
inverse scattering transform method for the DS equations was formulated 
in 0-0. This method allows to construct the soliton solutions and to 
study the global properties of the solutions, for instance, the asymptotic 
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behavior at large time [10, 11]. The asymptotic behavior of the solutions for 



nonintegrable DS equations was studied in [|T^ 



The perturbed equations (|T]) are nonintegrable by inverse scattering trans- 
form method. Explicit form of the perturbation operator can be defined by 
small irregularity of bottom or by taking into account the next corrections 
in more realistic models for the liquid surface considered in [jl], 0. As main 
example in this work we use the perturbation appearing because of small 
irregularity of the bottom: F = AQ. Here A is real constant. The sign (±) 
of A corresponds to decreasing or increasing depth with respect to spatial 
variable ^. Here we don't discuss the other perturbations corresponding to 
the system of liquid surface though the method we used allows to study the 
DS-1 equations with more wide class of the perturbation operators. 

The special solution of DS-1 we perturbed here was constructed in the 
work 0: 

^^ pXfiexpjitjX^ + fi'^)) 

^ 2cosh(//0cosh(Ar;)(l -a|p|2^(l + tanh(A7]))(l + tanh(/iO))' 

(2) 

where A, /x are positive constants defined by boundary conditions as ?] ^ — oo 
and ^ — oo; p is free complex parameter of the solution. 

This solution decreases with respect to spatial variables exponentially. 
Besides if cr = 1 and > 1 this solution has singularities on some lines 

rj = const and ^ = const. 

The solution was called dromion in work 0. The inverse scattering 
method for dromion-like solution of was developed in 0. 

The constructing of the the asymptotic soliton-like solutions for the per- 
turbations for the integrable equations is popular field. The perturbations 
of the (l-l-l)-dimensional integrable equations are investigated in more de- 
tail. In that case the asymptotic solutions are usable at long time and the 



Fourier- like integral corrections were studied in |I3|-|21]. The instability of 



the solutions of the (2-|-l)-dimensional equations with respect to transversal 



perturbations was studied for example in |^-|^. But the perturbation of 
the essential two-dimensional solutions with respect to spatial variables of the 
(2-1-1 )-dimensional integrable equations is less developed field. First cause is 
more diversity of the solutions and second cause is more bulky formulas for 



the solutions. Here we must refer the works [^-[0- In these works the 
perturbation of the soliton solution for the DS-2 equations had been studied. 
From these results and from the work about the asymptotic behavior of non- 
soliton solutions for the DS-2 equations it follows that the soliton of the 
DS-2 equations is instable with respect to small perturbation of the initial 
data. These results stimulate the studies of the dromion perturbation. 
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The possibility of the full investigation of the linearized DS-1 equations 
plays the main role to construct the perturbation theory of the nonlinear 
DS-1 equations. The linearizations of the integrable equations have basic set 
of solutions usually. It allows to solve the linearized equations by Fourier 
method. For (l+l)-dimensional equations it was shown by Kaup in [p^] . 
The works ||2^-[0 were devoted to the same questions. The basic sets of 
the solutions for the linearized DS-2 and DS-1 equations were obtained in 



1 Problem and result 

We construct the asymptotic solution of equations (|lD on mod(0(£^)) uni- 
formly at large t. The perturbation operator is F = AQ and the boundary 
conditions for Gi and G'2- 

^2 

^^'^-°^ = ^^"2cosh^(Ar^)' ^^'''-°° = ^^"2cosh^(/.0- 
We find the asymptotic solution in the form: 

Q(e, V, t, e) = W{^, T], t, t) + eU{C, V, t, r) , 
Gii^, V, t, e) = gi{C, 77, t, t) + ehi{^, 77, t, r), 
G2{^, V, ^) = 92{^, V, ^) + ^^2(?, V, (4) 

where t = et is slow time. The leading term of the asymptotics has the form: 

W{C,r],t,T) = q{^,r],t-p{T)), 

and Qi, g2 are: 

Z J —00 
z J —00 

Denote 7(r) = 1 - |p(t)|2, 70 = 7(0). 
Theorem 1. // 

7(^) = 70""^^^^^^ Arg{p{T)) = const, 

where > 1 at a = —1 and 0<'~fo<lata = l, then the asymptotic 
solution with respect to mod{0{e'^)) is useful uniformly at t = 0{e^^). 
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Remark 1. At more long time t ^ e"^ log(log(£:^^)) the formulas (Q) 
are asymptotic solution of (|l|) with respect to mod(o(l)) only. 

The asymptotic analysis given here is valid for the solutions without 
the singularities. It means if a = 1, then ^|p(t)P < 1 or the same 7(r) > 0. 
If the coefficient of the perturbation A > 0, then \p\ increases (7 — > at r ^ 
00) with respect to slow time. It allows to say that the singularity may appear 
in the leading term of the asymptotics at r ^ 00. However we can't state 
this strongly because our asymptotics is usable only at r -C log(log(£~^)), 
while the order of j{t) is greater than e. The order of a remainder term after 
substituting the asymptotic solution into the equations (P is only o(l) 
at the very long time 1 ^ r ^ log(log(e~^)). 

The obtained result for the problem about the interaction of the long and 
short waves on the liquid surface shows that when the depth decreases the 
the formal asymptotic solution can be described by adiabatic perturbation 
theory of the dromion for the focusing DS-1 equations at least for \t\ -C 

£-Mog(log(£-l)). 

Remark 2. The appearance of the singularities in the solution of non- 
integrable cases of the Davey-Stewartson equations is known phenomenon 



33. 



Remark 3. D. Pelinovsky notes to author, that the modulation of the 



parameter \p\ may be obtained out of the "energetic equality" 
dt I L didr^\Q\^ = e [ L d^dr]{QF - QF). 



2 Solution of linearized equations 

Here we remind the formulas for solution of the linearized DS-1 equations on 
the dromion background: 

idtU + (a| + dl)U + (Gi + G2)U+ (Gi + G2)Q = iF (5) 
a^Fi = -|a,(Qt/ + QU), d,V2 = ^d^iQU + QU). 

The results of inverse scattering transform O and the set of the basic 



functions |33| are used for solving of the linearized equations by Fourier 



method. However unlike the work here the solution of the DS-1 equations 
with nonzero boundary conditions is considered. It leads to the change of 
the dependency of the scattering data with respect to time (see also 0) and 
of the formulas which define the dependency of Fourier coefficients of the 



solution for the linearized DS-1 equations in contrast to obtained in |]33 |. 
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In the inverse scattering transform one use the matrix solution of the 
Dirac system to solve the DS-1 equation 0-0: 

Let and i/j' are the matrix solutions of the Goursat problem for the 
Dirac system with the boundary conditions (M) 



(7) 



Denote by -0^-), j = 1, 2, the column of the matrix then this column is 
the solution of two systems. There are the system (|) and additional system 
of time evolution: 



One can obtain the equation like this for the other columns of the matrices 
These equations are differ from (j^) by the sign of ik"^ in first term of the 
right hand side only. 

Now we write two bilinear forms defining analogs of the direct and inverse 
Fourier transforms. First bilinear form is 

/OO /"OO _ 
/ dCdr]{xifii(rf + X2fJ'2f)- (9) 

Here Xi ^-^^d fii are the elements of the columns x ^-^^d /i. 

Denote by and (/)(^2) the solutions conjugated to ^^-^ and 'ifj'^2) with 
respect to the bilinear form (|^). 

Using the formulas for the scattering data ([0]) one can write these data: 

s,{k,l) = i-(^+)(e,r7,A;),E(i)(z;0)Q, (10) 

^2(^,0 = ^(^(2)(^'^'^)'^(2)M)q- (11) 

Here E{z) = diag(exp(2;), exp(— 2;)). 
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It is shown ([0]) the elements of the matrices are analytic functions 
with respect to the variable k as ±Im{k) > 0. Using the scattering data one 
can write the nonlocal Riemann-Hilbert problem for the ipii and 1^12 on the 
real axes (0): 



i'lii^^V^k) = exp{ikr]) + exp{ikr]) i^exp{-ikr]) J^^dl Si{k,l)^pf2i^,r],l) 
ipi2{^,V,k) = exp{-ikC)(^exp{ik^) J^^dls2{k,l)'i/j^^{C,r],l)^ . 

Here 



2i7T J -00 k' -{k± iO) 
The Riemann-Hilbert problem for -0^^ and -022 has the form: 



^2i(?>^5^) = exp(iA;7])(^exp(-z/i;r/) /^d/ Si (A;, 0^22(^,^,0 
V'22(?,^,^)= exp(-2A;0 + 

exp(-?A;^) expiik^) dl S2{k, 0V'2i(?, V, 

Introduce second bilinear form 

{X, l^)s= / dkdl{x\l)f^\k) S2{k, I) + x'{l)f^'{k) s,{k, I)), (12) 



00 J —QO 



where X"' is the element of the row x- 

Denote by ip^^\ j = 1,2, the row conjugated to ip^^^ = [tpji,ipj2] with 
respect to the bilinear form (12). Formulate the result about the decompo- 
sition obtained in |]53[ . 

Theorem 2. Let Q be such that d"Q e LiDC^ for \a\ < 3, if f{^,v) is 
(9"/ G Li n for \a\ < A, then one can write f in the form 

f = —{i^^'\^,V,l),V^'\^,V,k))f, 
where ^ 

Using the theorem 2 one can solve the Cauchy problem for the linearized 
DS-1 equations. Here the dependence of / with respect to t is differ from the 
same obtained in because here we use the DS-1 equations with nontrivial 
boundary conditions (|). These changes may be obtained using the results 
of § and 
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Theorem 3. Let Q be the solution of the DS-1 equations with the boundary 
conditions Gi|g^_oo = and G2\r^^-oo = ^2, and Q satisfy the conditions of 
theorem 2, the solution of the first of the linearized DS-1 equation is smooth 
and integrable function U with respect to ^ and rj, where d°'U G Li fl and 
e Li n C\ for \a\ < A and t e [0,To]. Then 

POO 

dtU = i{k^ + l^)U + / dk'U{k - k', I, t)x{k') + 

dl'Uik,l-l',t)K{l') + F, (13) 

-oo 



If the boundary conditions in the problem for the DS-1 equation are zero, 
then X = 1^ = 0- In this case the formulas of the theorem 3 allow to solve the 
linearized DS-1 equation in the explicit form. In this work we consider the 
solution of the DS-1 equation with nonzero boundary conditions. It leads to 
integral terms in the formulas of theorem 3. In order to solve the linearized 
DS-1 equation we must transform the formula (0). In the right hand side 
of (p!3D the integral terms are the convolutions. Go over to the equations for 
the Fourier transform of the functions U{k, l,t,T) with respect to variables k 
and I. As result we obtain the linear Schrodinger equation: 



Here 



idtU + {dl + dl)U + {U2{i^i) + ui{ri\))U = F. (14) 

U (^, ri,t,T) = — [ dkdlU{k, I, t, r) exp{—ikri — il^); 

F{^, ri,t,T) = — / , dkdlF{k, I, t, r) exp{—ikrj — ilC,)- 
2% JR^ 

The same equation without the right hand side was obtained in [0] for the 
time evolution of the scattering data for the nonlinear DS-1 equation. 

One can obtain the solution of the Cauchy problem with zero in the initial 



conditions for the equations (|T^ by the separation of the variables. In our 
case the solution of the equations (|T^ obtained by the Fourier method has 
the form: 

U{^,ri,t) = — f drndnij (m, n, t)X{^, m)Y{r], n) exp(— zt(m^ -|- n^)) + 



2ti 



1 



V27r 
1 



dnU^{n, t)Y{7], n)X^{^) exp{-it{n^ - /x^)) + 
dntlxim, t)X{^, m)Yx{r]) exp{-it{m^ - A^)) + 
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Here we use notations: 

X{m, = — : exp(-zmO, ^^(6 



im — jj ' 2cosh(/z^) 

Atanh(A?7) + in r ■ ^ \ ^ 

Y{n,r]) = : exp{-mr]), Yx{r]) - 



in — X ' 2cosh(A?]) ' 

dtU = i[m^ + n^)?/ + F(m, n, t), dt% = i{r? - f/)U + F^{n, t), 
dtlJx = %{jr? - X')U + h{m, t), dtU^x = -t{X^ - i?)U + 

U\t=o = Ufj,\t=o = U\\t=o = Ufj,\\t=o = 0; 



F(m, n, t) = d^dr]F{^, r], t)X{m, C)Y{n, 7]), 

Fxim,t,T) = j didr,F{i,'n,t)X{m,i)YM, 
F,x{t,T) = j^Jidr^F{i,r^,t)X^{OYM- 

3 The equation for first correction 

In this part the equation for the slow modulation of the parameter p(r) is 
obtained. This equation is necessary and sufficient condition for the uniform 
boundedness of the first correction of the expansion (g). 

Substitute the formula (^) into the equations (|l]). Equate the coefficients 
with the same power of e. The equations at are realized because W,gi,g2 
are the asymptotic solution of nonperturbed DS-1 equations. For the first 
correction we obtain the linearized DS-1 equations: 

idtU + (9| + d^)U + {gi + g2)U + {Vi + V2)W = iH 
d^Vi = -^drtiWU + WU), d^V2 = ^d^{WU + WU), (15) 

where 

H = AW- drW. 

Before to use the formulas from the previously section we reduce the form 
of the right hand side in first of the equations ([T5|). In the leading term the 
parameter p depends on the r only. The other parameters depend only on the 
boundary conditions and do not change under perturbation. For convenience 



8 



we represent p(r) = r(T) exp(ia(r)). The derivation of W with respect to 
slow variable r can be written as: 

d^W = drWr' + d^Wa'. 

Here the derivatives r' and a' are unknown. They will be obtain below. 

Compute the function H. From the theorem 2 and formulas of the func- 
tions -0+ and (j) (see Appendix) we obtain: 

H{k, /, t, r) = exp(-it(A2 + /x^)) {p{k, I; p) - R{k, I; p)drp^ , 

where 

P{k, I; p) = exp{it{X'^ + p'^))AW, R{k, I; p) = exp(it(A2 + p^))drW. 

In these formulas we write the dependence on time in explicit form. It allows 
to remove the secular terms in the asymptotic solution (^. The differential 
equations for U have the forms: 

dtU = i{m^ + n^)U + exp{—it{X^ + p?)){P{m, n; p) — R{m, n; p)), 
dttj^ = t{n^ - p^)U + exp(-2t(A2 + p^)){P^{n; p) - R^{n; p)), 
dtUx = i{m^ - X^)U + exp(-^t(A' + p^))iPxim; p) - Rx{m- p)), 
dtU,x = -^{X^ - f^')U + expi-ttiX' + p^)){P,x{p) - RM)- 

One can see that the secular terms may appear because of the last term in 
the equation for U^x ■ The requirement of equivalent to zero of this term leads 
us to the equation for p(r): 

-R^A - = 0, p\t=o = po- (16) 

In the result Ufj_x = 0. The other equations for (j are easy to integrate. 
The solutions of this equations are bounded on the all of arguments for all 
time. 

We must return to the original of the images U to state about bounded- 
ness of the solution U, Vi, V2 for the equations (^). One can see the direct 
(from U into U) and inverse (from U into U) integral transforms as the 
Fourier transform from the smooth and exponentially decreasing functions 
with respect to the corresponding variables. The Fourier transform moves 
such functions into analytic functions near the real axis. The inverse trans- 
form moves these analytic functions into the exponential decreasing func- 
tions. So the solution of (|15|) is boundedness and decreasing exponentially 
with respect to the spatial variables. 
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4 Modulation equation for p(r) 

Here the equation ([T6|) for the parameter p{t) is reduced to the more convi- 
ent form. Write the derivative of the leading term with respect to the slow 
variable r. 

r' 2iW r' 

id^W = -a'W - iW- + 



r ' 1 - (Tr2^(l +tanh(^0)(l + tanh(A?7)) r 
Denote F = and compute the images {■)^x every term. 

(W)^, = -Ji:^2tfl±I^^,r - 1) log 11 - .ri; 

Denote the image of the last term by /i^a- Its image has the form: 

V = ^ (1 - -r) (^^-^ - 1 - log |1 - ari 

The image {■) of AW has the similar form. 

Substitute these formulas into ([T6| ) and separate real and imaginary parts 
of this equation, then 

a' = 0, 

-{aV - 1) log |1 - ar| - - far - (1 - aV) log |1 - aT\) + 
r r \ J 

A(ar-l)log|l-ar| = 0. 
Use the notation for F, then the second equation has the form: 

— = -2aA{l - aT) log |1 - ar\. 
dr 

This equation defines the evolution of the absolute value of the complex 
parameter p. The argument of this parameter do not change under the per- 
turbation F = AQ. 

Solve the equation for F. Denote 7=1 — aT, rewrite the equation for 7, 
then we obtain: 

7' = 2A7log(7). 
The solution for this equation has the form: 

7(r) = exp(Cexp(2Ar)), 
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where 7|r=o = Gxp(C), then we can write the 7(r) in the form: 
The theorem 3 is proved. 
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5 Appendix 

5.1 Explicit formulas 

Here we remain the explicit forms of the solution for the Dirac equation 
with the dromion-like potential. These forms were obtained in 0. In our 
computations we use first column of the matrix -0+ only. 

\ _ exp(zfcr^) \ Jri 2cosh{Ap) 



V ^^21 ; V / (l-a|p|2^(l + tanh(/iO)(l + tanh(A77))) 

/ -a|ppAM(l+tanh(pg)) \ 

8cosh{An) /i-vN 
^ I apMexp(-»f(A2+p^)) • y^') 
\ 2cosh(M0 / 

5.2 Solutions of conjugated equations 

Here we write the problems for the functions conjugated to -0^ with respect 
to the bilinear forms. 

The functions are the solutions of the boundary problem conjugated to 
the solutions of the problem with respect to the bilinear form (|^). 

First column of the matrix is the solution of the integral equation: 

Z J —oo 

—1 r°° 

<P2i{Lv,l,t) = — dCQ{e,v,t)Me,v,l,t). 

The explicit formulas for first column of the matrix used in section 4 
for the dromion potential has the form: 

( '^ll 1 = ( ^^P(^^^) 1 I -ig 2cosh(pp) ^ 

; V / (1 - ^IpPtKI + tanh(A.O)(l + tanh(Ar;))) 
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-g-|ppA^(l+tanh(AT?)) 
/ 8cosh(^tg) 
^ ' -o-pAexp(-»f(A^+/j^)) 
2 cosh(A»7) 



Second bilinear form (jT^) allows to write the integral equations conju- 
gated to integral equations which were obtained from the nonlocal Riemann- 
Hilbert equation for the functions tp"^ in 0. Using this equations for the 
functions r], I) one can show that the functions are the solutions of the 
boundary equations for the Dirac system or the integral equations: 

1 p 

Z J —oo 
L J —oo 



The functions ^p have the similar form as the functions 
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